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Abstract 

The electromagnetic pion production reactions are investigated within the dynamical coupled- 
channels model developed in Physics Reports, 439, 193 (2007). The meson-baryon channels 
included in this study are jN, ttN, r]N, and the irA, pN and aN resonant components of the 
irnN channel. With the hadronic parameters of the model determined in a recent study of ttN 
scattering, we show that the pion photoproduction data up to the second resonance region can be 
described to a very large extent by only adjusting the bare 7-/V — * N* helicity amplitudes, while 
the non-resonant electromagnetic couplings are taken from previous works. It is found that the 
coupled-channels effects can contribute about 30 - 40 % of the production cross sections in the A 
(1232) resonance region, and can drastically change the magnitude and shape of the cross sections 
in the second resonance region. The importance of the off-shell effects in a dynamical approach is 
also demonstrated. The meson cloud effects as well as the coupled-channels contributions to the 
7./V — ► N* form factors are found to be mainly in the low Q 2 region. For the magnetic Ml — > A 
(1232) form factor, the results are close to that of the Sato-Lee Model. Necessary improvements 
to the model and future developments are discussed. 
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I. INTRODUCTION 



It is well recognized [j], 0] that the data of electromagnetic meson production reactions 
can be used to reveal the structure of the excited states (N*) of the nucleon. In this paper, 
we report on an investigation in this direction within the dynamical coupled-channels model 
(MSL) presented in Ref. [3| which is being applied at the Excited Baryon Analysis Center 
(EBAC) of Jefferson Laboratory. 

The coupled- channels approach has been used 0, 0, @, 0] in recent years to analyze the 
meson production reaction data. It is therefore useful to briefly emphasize here the essence 
of taking a dy namical _approach as developed in Ref. jl] and in earlier works @, 0, 10, 11, H, 
H, HQ EH HQ. Since N* states are unstable, their structure must couple with the 
reaction channels in the meson production reactions. To determine correctly the spectrum of 
N* states, an analysis of the meson production data must account for the coupled-channels 
unitary condition. The extracted N* parameters can be interpreted correctly only when the 
reaction mechanisms in the short-range region, where we want to map out the N* structure, 
have been accounted for. The MSL model meets these two crucial requirements and is 
therefore suitable for analyzing the world data of meson production reactions induced by 
pions, photons, and electrons. 

Schematically, the MSL model solves the following coupled integral equations in each 
partial wave 

T a , P (Pa,P(3; E) = V aj p(p a ,Pf3) +J2 P 2d P V a,5(Pa,p)G S (p, E)T S;f3 (p, pp, E) , (1) 

<5 J 

with 

V a ,p(p a ,Pl3) = V a ,p{p a ,pp) + 2^ E-M% ' ( ' 

where a,(3,5 = jN, ttN, rjN, and 7T7riV which has 7rA, pN, aN resonant components, G$(p, E) 
is the propagator of channel S, M^* is the mass of a bare excited nucleon state iV*, v a> p is 
defined by meson-exchange mechanisms, and the iV* — > (5 vertex interaction T^*,/? is related 
to the quark-gluon sub-structure of iV*. If we take the on-shell approximation, Eq. (pQ) is 
reduced to the following algebraic form of K-matrix models 



T lfi(P^Pfi, E ) = V aAPa>Ps) x [85,0 + ip(Ps)Ts^(p s ,Pf3, E)} , (3) 

8 

where p(ps) is an appropriate phase space factor. Qualitatively speaking, models that make 
use of on-shell expressions such as Eq. fl3]) are used to avoid an explicit treatment of the 
reaction mechanisms in the short range region where we want to map out the quark-gluon 
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FIG. 1: The non-resonant electromagnetic meson production interaction vmb^n, where MB 
7rN,r]N,TrA,pN,aN. The details are given in Ref. [3(. 
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sub-structure of N* states. Thus the N* parameters extracted by using Eq. (pQ) can be more 
directly interpreted in terms of the quark-gluon sub-structure of iV*. From the study [ll, ill 
in the A (1232) region, it is reasonable to interpret r^v*,/? in terms of hadron structure 
calculations with effective degrees of freedom, such as the constituent quark model [23J and 
the model [lij based on Dyson-Schwinger Equations. In the near future, one hopes to relate 
Fn*,/3 to Lattice QCD (LQCD) calculations |25l . 26[. The possibility of making contact with 
the current hadron structure calculations, which can be carried out with sufficient accuracy 
in the forseeable future, is the main motivation for performing the analysis using a technically 
much more involved dynamical approach. 

To analyze the electromagnetic meson production reactions within the MSL model, it is 
necessary to first determine the hadronic part of its Hamiltonian. A progress in this direction 
has been made recently in Ref. [27| (JLMS). The main purpose of this work is to explore the 
consequence of the parameters determined by JLMS in describing the pion photoproduction 
reactions. Within the MSL formulation, these hadronic parameters should be consistently 
used to fix the strong interaction vertices of the non- resonant 7 A — > MB interaction vmb^n 
of Eq. (T5]). Since the electromagnetic interactions, such as the •yNN^np, •ynuj vertices of 
the considered non-resonant mechanisms illustrated in Fig. HJ have been determined in the 
previous works as also given in Ref. 0], the challenge here is to explore whether the pion 
photoproduction data can be described by only adjusting the bare — > N* vertex Tn'^n 
ofEq. ©. 

To proceed, we first note that in the study by JLMS j27[ it was found that the fit to ttN 
elastic scattering data is not sufficient for pinning down precisely the hadronic parameters 
associated with the unstable particle channels 7rA, pN, and aN . These channels have 
very large effects in the invariant mass W > about 1.65 GeV region where the two-pion 
production dominates. It is necessary to include the irN —>■ tckN data in the fit. This 
very challenging task is still being pursued at EBAC For this reason, we will limit our 
investigation to the W — 1.1 GeV — 1.65 GeV region which covers the energies of the low- 
lying nucleon resonances in the so called first and second resonance regions. The resulting 
model is sufficient for investigating the dynamical coupled-channels effects on the 7jV — > nN 
cross sections and the jN — > N* transitions. This is the main focus of this work. 

The details of the employed dynamical coupled-channels model have been given in Ref. 0] . 
In section II we only recall the formulas relevant to the pion photoproduction reactions. 
The results are presented in section III. Section IV is devoted to discussing the necessary 
improvements and future developments. 



II. FORMULATION 

In the helicity-LSJ mixed-representation [§] where the initial 7A state is specified by its 
helicities A 7 and A at and the final MB states by the (LS)J angular momentum variables, the 
reaction amplitude of 7(g) + N(—q) — > 7i(k) + N(—k) at invariant mass W can be written 
within the MSL formulation as (suppress the isospin quantum numbers) 

(k,q,W)=t 

L Sj\[ 7T TV , A -y A jv 

(k,q,W)+t (k,q,W), (4) 

where S^- = 1/2 is the nucleon spin, and the non-resonant amplitude is 

^lsnttN^Xn^i^i^) = v ls n ttN,x 1 x n (^^ 9) E) + ^2 ^2 k' 2 dk't'L SN7TNLIS , M/BI (k,k f ,E) 

M'B' US' 
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xG MIBI (k',E)vl, s , MIB ,^ XN (k',q,E) . (5) 

In the above equation, the meson-baryon channels included in the sum are M'B' = 
Ti N, rjN, ii A, pN, o~N. The matrix elements v[ SMBX \ N {k,q, E), which describe the — > 
M'B' transitions, are calculated from the tree-diagrams, illustrated in Fig. [TJ of a set of 
phenomenological Lagrangians describing the interactions between 7, it, rj, p, u, a, N, and 
A(1232) fields. The details are given explictly in Appendix F of Ref. [3(. We will use the 
hadronic parameters determined by JLMS |27| to evaluate the meson-baryon propagators 
Gm'b' and the nN — > MB amplitudes t"[ s N L , s , M , B ,(k, k', E), and to also fix the hadronic 
vertices of the — > MB amplitudes v[ SMBX ^ XN (k,q,E). As discussed in section I, all 
of the electromagnetic vertices, such as jNN, "fmr, ■yTru, in v[ SMBX x (k,q,E) are taken 
from previous works, as also specified in Ref. [3j. Thus the non-resonant amplitude defined 
by Eq. (jSJ) is completely fixed in the present investigation. Such a consistent dynamical 
treatment of strong and electromagnetic reaction mechanisms is highly desirable in using 
the meson production reactions to study the N* structure. 
The resonant amplitude in Eq. is 

tL's N ,N,x,x N (k,q,E)= £ K w (^)rA,iWfi., v J#), (6) 

n;,n; 

where the dressed vertex functions are defined as 



[^N*,LS N irN{k, W)]* = [T J N ^ LSNnN (k)}* 



(k, k', W)G M , B ,{k\ W)Kwm>b>W)\* 

H IT I Tit T I t~tt *J 

(7) 



M'B' L'S' 



rAr* ) A 7 A J v(?) 



+ E E / k' 2 dk'T 

N*, L'S' M'B' 

(k f , W)G MIB ,(k', W)v J L , s , MIB ,^ XN {k', q) .(8) 



M'B' L'S 



The second term of Eq. (jSJ) is due to the mechanism where the non-resonant electromagnetic 
meson production takes place before the dressed N* states are formed. The contribution 
due to the ttN intermediate state is illustrated in Fig. [2j Similar to what was defined in 
Ref. 0,13, we call this contribution the meson cloud effect to define precisely what will 
be presented in section III. 

The N* propagator in Eq. ([6]) is defined by 



[D(W)-% = (W- M°,)5 M - , 
where M^* is the bare mass of the N* state, and the self-energies are 

(k)G MB (k,W) 



MB LS 



dkk 2 Y J Nf L SM L > \ 



1 N*,LSMB 

(k,W) 



(9) 



(10) 



Since the dressed vertex T J N * LSMB {k, W) of A^* -> MB in Eqs. ©-© and the bare 
masses Mjy* and the mass shifts Ejj in Eq. ([9} have been determined by JLMS [27|, the 
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FIG. 2: Graphical illustration of the contribution to the ttN intermediate state to the dressed 
7iV N* vertex defined by Eq. 

only unknown quantities in the resonant amplitude Eq. ([6]) are the bare 7iV — > N* vertex 
functions r^*^ \ N (q) in Eq. ([8]). We parameterize these functions as 



rj.^fa) - ^^y^^T7^[\/2^ A A]^(A T -A N ), (ii) 

where g# is defined by the N* mass M^* = qn + Ejsr(qii)- For later discussions, we also cast 
the dressed vertex Eq. ([8]) into the form of Eq. (jTTJ) with a dressed helicity amplitude 

A{{q)=A{ + AT^q) (12) 

where A™ C,J \q) is due to the meson cloud effects defined by the second term of Eq. ([H]). The 
dressed helicity amplitude A J X is related to the partial width of the electromagnetic decay of 
a N* resonance defined by 

[partial decay width] = f|^^^p3/ 2 M| 2 + \A J 1/2 (q R )\ 2 ] . (13) 

Eqs. (TTTl)- (iT3l and Eq. fix the relation between the bare helicity amplitude A J X of Eq. ([8]) 
and the helicity amplitudes A J X listed by the Particle Data Group (PDG). 

III. RESULTS 

With the specifications given in section II, our task is to perform x 2 — fits to the available 
data of 7iV — > ttN reactions up to W = 1.65 GeV by adjusting the bare helicity amplitudes 
A x of the bare iV* states included in the JLMS fit [27] to the ttN scattering data. The 
X 2 -fits are performed by using MINUIT. In this first and second resonance regions, we find 
that the fits to the data are mainly sensitive to the bare helicity amplitudes listed in Table 
HI Other bare helicity amplitudes are set to zero in the calculations. The quality of the 
resulting fit can be seen in Figs. [31 SI El and The solid curves in Figs. [3] and H] are the 
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TABLE I: The bare 7 A — > A* helicity amplitudes determined from x 2 -fits to the 7 A — > ttN data 
shown in Figs. [3][6l 



calculated differential cross sections (da/dfl ) for 7£> — > 7r°£> and 7^ — > 7r + n, respectively, 
compared with the corresponding experimental data. Similarly solid curves in Figs. [5] and 
[6] correspond to the obtained photon asymmetry (E 7 ) compared to the experimental data 
for the reactions 7^ — > 7r°p and 7^ — > 7r + n, respectively. We see that the model can give 
an overall good description of the considered data, while significant discrepancies with the 
data remain. 

We emphasize here that the determined bare helicity amplitudes listed in Table |T] are 
not directly the properties associated with the nucleon resonances. They are simply the 
properties of the excited nucleon states in the absence of coupling to the reaction channels. 
We need to identify the resonant positions from the partial-wave amplitudes predicted by our 
model. The dressed helicity amplitudes Aj^ calculated at those resonance positions according 
to Eq. (jSJ) can then be compared to the partial decay widths from the measurements, as 
seen in Eq. ffl~3]) . This is being pursued by developing j29| an analytic continuation method 
to evaluate the reaction amplitudes in the complex energy plane following the dynamical 
coupled-channels equations of the MSL model [3J. 

We now turn to investigating the coupled-channels effects. In Fig. we see that the 
calculated total cross sections (solid curves) are in good agreement with the data. The 
dashed curves are obtained when the channels nN, it A, pN, and aN are turned off in the 
loop integrations in the non-resonant amplitude defined by Eq. (jSJ) and the dressed 7^ — > A^* 
vertex defined by Eq. ([8]). Clearly, the coupled-channels effects 7^ — > nN, it A, pN, aN 
— > 71N can change the cross sections by about 10 - 20 % in the A (1232) region and as much 
as 50 % in the W >1400 MeV second resonance region. The corresponding coupled-channels 
effects on the differential cross sections are shown in Fig. [HI Comparing the solid and dashed 
curves, we see that the coupled-channels effects can change the magnitudes and shapes of 
da/dfl, in particular at higher energies. 

As discussed in the introduction, an essential feature of a dynamical approach is to solve 
integral equations Eqs. ([I])-® which involve off-shell amplitudes. The off-shell effects are 
due to the reaction mechanisms at short distances which are not treated explicitly in the 
technically much simpler K-matrix coupled- channels models 0, 0, @, 0,0,0. If the off- 
shell effects are neglected in the 7rA^-loop only calculations, we then bring the dashed curves 
in FigsJT] and [8] to the dotted curves. Clearly, the off-shell effects are very significant, as was 
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FIG. 3: Differential cross section for 7p — ► ir°p compared to experimental data obtained from 
Ref. Q. 



also revealed in the coupled- channels calculations [l8[ of KY photoproduction. 

To disentangle the structure of nucleon resonances, it is important to investigate the 
meson cloud effect on the 'jN — > N*, as defined by the second term of Eq. (jSj). It can have 
contributions from the loop integrations over the ttN, tjN, it A, pN, aN intermediate states 
in this calculation. The contribution from the ttN loop is that illustrated in Fig. [2l If the 
electromagnetic form factors associated with the non-resonant interaction vmb^n are taken 
from previous works, we can predict the Q 2 -d e P en dence G f the me son cloud effect term A™ C ' J 
ofEq. (P). 

We first investigate the 'jN — > A (1232) transition. The resonance position of A (1232) is 
well defined and reproduced in our calculations as can be seen in the good description of the 
cross section data near W = 1232 MeV in Figs. [3H3 For this isolated resonance, we can follow 
the procedures detailed in Ref. llj to calculate the 'jN — > A (1232) magnetic form factor 



G* M {Q 2 ) from the imaginary part of the multipole amplitude M ± + of 7*iV — > ttN reactions 
at W = 1232 MeV. Our results are shown in Fig. O The dashed curve corresponds to the 
full meson cloud contribution from this calculation. The dotted curve is obtained by keeping 
only the 7r-loop in the calculation. Clearly the difference between the dashed and dotted 
curves is due to the coupled-channels effects 7iV — > i]N,irA, pN,aN — ► A (1232). Here we 
note that the pion-loop only result (dotted curve) is very close to the meson cloud effect 
predicted by the Sato-Lee Modelpj]]. The solid curve following the data is just for guiding 
the eyes. The difference between the solid curve and the dashed curve provide information 
about the Q 2 -dependence of the bare 7^ — > A (1232) form factor which can be used as the 
starting point of our dynamical coupled-channels analysis of pion electroproduction. 

For the meson cloud effects on the other 7^ — > N* vertex, we calculate the second term 
of Eq. © to get A™ C > J by using the normalization defined by Eqs. f|TTl) - ffT2l . In Fig. [TU1 we 



show the predicted mag nitudes of \A™ C ' J \ evaluated at W = 1535 MeV for S n , W = 1440 
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FIG. 4: Differential cross section for — > 7r + n compared to experimental data obtained from 
Ref. 0. 



MeV for P u , W = 1520 MeV for D 13 , W = 1625 MeV for D 15 , and W = 1620 MeV for 
. The solid dots at Q 2 = are the determined bare helicity amplitudes. The solid curves 
are from the full calculations and the dashed curves are from only keeping the ttN loop in 
Eq. (jHJ). We see that the meson cloud contributions and the coupled-channels effects affect 
mainly the low Q 2 region. 

Here we note that the results presented in Fig. [TD] are around the resonance positions 
listed by PDG, not from the resonance pole positions which will be determined 2^] using an 
analytical continuation 29]. Thus the results presented here are only for giving some qualita- 
tive estimate of the meson cloud effects on 7./V — > N* excitation. More accurate predictions 
will be published in our subsequent analysis [1^] of the data of pion electroproduction. 



IV. SUMMARY AND FUTURE DEVELOPMENTS 



We have applied the dynamical coupled- channels model of Ref. [3J to investigate the 
pion photoproduction reactions in the first and second nucleon resonance region. With the 
hadronic parameters of the model determined in a recent study [27] of nN scattering data 
and the non-resonant electromagnetic couplings taken from the previous works, we show 
that the available data of differential cross sections and photon asymmetries of 7V — > ttN 
up to W = 1.65 GeV can be described to a very large extent. The only free parameters in the 
X 2 -fit to the data are the bare 7V — > N* helicity amplitudes. It is found that the coupled- 
channels effects can have about 10 - 20 % effects in the A (1232) resonance region, and can 
drastically change the magnitudes and shapes of the cross sections in the second resonance 
region. We also demonstrate the importance of the off-shell effects in a dynamical approach. 
The meson cloud contributions to the jN —>■ N* form factors have been predicted. For all 
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cases, they are mainly in the low Q 2 region. For the magnetic Ml 'yN — > A (1232) form 
factor, the results are close to those of the Sato-Lee Model. The coupled-channels effects on 
the meson cloud contributions are also found to be mainly in the low Q 2 region. 

The results presented here can be the starting point for analyzing the pion electropro- 
duction data in the first and second resonance region. It will be interesting to see the bare 
7iV — > A (1232) form factor which can be estimated from the differences between the solid 
and dashed curves in Fig. [9] can be verified in the analysis of the data at W = 1232 MeV 
and Q 2 < 6 (GeV/c) 2 . For the data in the second resonance region, our task is to extract 
the bare helicity amplitudes of 7 iV — > N* at each Q 2 , as was done in this work at the photon 
point Q 2 = 0. Of course the data to be fitted are much more extensive and we need to also 
determine the longitudinal components of the jN — > iV* vertex. This is being pursued [3(| 
at EBAC. 

The most unsatisfactory part of this work is the uncertainties in determining the bare 
7iV — > N* helicity amplitudes. As seen in Table [U there are 16 helicity amplitudes adjusted 
in our x 2_ fits to the data. Obviously, it is possible that there could exist other solutions 
to the minimization problem if more experimental data, in particular polarization data is 
included in the fit. Within the MSL formulation 0], these bare parameters are related to 
hadron structure calculations in the absence of the couplings with reaction channels. A more 
fruitful approach is to take the helicity amplitudes predicted from such hadron structure 
calculations as the starting values in our x 2_ fit to avoid unphysical solutions. The resulting 
parameters could then be used to examine that hadron structure calculation. However 
more theoretical analysis is needed to know precisely what kind of structure calculations are 
consistent with the MSL formulation and can be used for this purpose. Since it was found 1^ 



that the extracted bare 'jN — > A (1232) magnetic Ml form factor is fairly consistent with the 
prediction of constituent quark models, one possibility is to use the relativistic constituent 
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FIG. 7: Total cross sections. The dashed curves are obtained from turning off all MB channels 
except the irN channel in the the loop integrations in the non-resonant amplitude defined by Eq. © 
and the dressed ^yN — * N* vertex defined by Eq. ([5]). The dotted curve is obtained by neglecting 
the off shell effects in the ttN only calculation. Experimental data obtained using Ref. 28l |. The 
diamonds correspond to the SAID solution Ec* 



quark model [31 



To improve the agreement with the data, one necessary next step is to improve the 
hadronic parameters of the model. These parameters, fixed at the values from the fit to 
ttN scattering data, must be improved by performing a combined analysis of both the ttN 
elastic scattering and ttN — > ttitN data. We expect that the parameters associated with the 
unstable particle channels, 7rA, pN, and aN, will be refined most significantly. Thus the 
predicted coupled-channels effects on 'yN — > ttN cross sections, as shown in Fig. [8] will be 
changed such that the fits to the data shown in Figs. [3HS1 in particular in the high W region, 
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FIG. 8: coupled-channels effects on da/dCl at W = 1.23, 1.30, 1.54 GeV. The dashed curves are 
obtained from turning off all MB channels except the irN channel in the the loop integrations 
in the non-resonant amplitude defined by Eq. ([5]) and the dressed 7iV~ — > N* vertex defined by 
Eq. ([8]). The dotted curve is obtained by neglecting the off shell effects in the 7riV only calculation. 
Right-hand side — ► 7r°p, left-hand side :jp — ► 7r + n 



can be improved. Furthermore, this is also needed to extend our investigation to the third 
resonance region where the two pion production dominates and the coupled-channels effects 
through these unstable particle channels are expected to be very large. 

Finally, we would like to address the questions concerning how the results from JLMS and 
this investigation can be used to extract the positions and widths of nucleon resonances. The 
resonances positions are the poles of the reaction amplitudes in the complex energy plane. 
The residues of these poles can then be related to the partial decay widths. If these poles are 
identical to the zeros of the N* propagator Dij(E) Eq. (jHJ), we then have the most desirable 
interpretation that the nucleon resonance is due to the coupling of the bare N* with the 
reaction channels. The dressed vertex functions Tn*,- 1 n, defined by Eq. (jSj), evaluated at 
these poles can be used to predict the dressed helicity amplitudes A\ for calculating the 
partial decay widths using Eq. ffl~3l) . On the other hand, the poles could be from the non- 
resonant amplitude such as the term t J LSN7rNX ^ XN of Eq. Then the identified resonances 
have nothing to do with the bare N* states and are simply due to the attractive meson- 
baryon interactions. Extraction of this resonance information requires developing numerical 
methods for solving the dynamical coupled-channels equations on the complex energy plane. 



Our effort in this direction will be published elsewhere [29 
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FIG. 9: The predicted Q 2 -dependence of meson cloud contributions to the magnetic Ml form factor 
G* M defined in Ref. Q and G D {Q 2 ) = 1/[1 + Q 2 /(0.71{GeV/c) 2 )} 2 . The dashed curve in the lower 
part of the figure is the predicted full meson cloud contribution. The dotted curve is obtained 
after turning off the coupled-channels effects due to irA and pN channels in the loop integrations 
of Eq. ([5]) . The data are compiled in Ref. [l9| . 



Acknowledgments 

We would like to thank Dick Arndt for useful discussions. We thank Mark W. Paris for 
pointing out a coding error in the calculation of the 7A— > pN matrix element. B.J-D 
acknowledges the support of the Japanese Society for the Promotion of Science (JSPS), 
grant number: PE 07021. B.J-D. thanks the nuclear theory group at Osaka University for 
their warm hospitality. This work is supported by the U.S. Department of Energy, Office 
of Nuclear Physics Division, under contract No. DE-AC02-06CH11357, and Contract No. 
DE-AC05-060R23177 under which Jefferson Science Associates operates Jefferson Lab, and 
by the Japan Society for the Promotion of Science, Grant-in-Aid for Scientific Research(c) 
15540275. This work is also partially supported by Grant No. FIS2005-03142 from MEC 
(Spain) and FEDER and European Hadron Physics Project RII3-CT-2004-506078. 



12 



200 




2 4 2 4 

S 2 [(GeV/c) 1 ] Q 1 [(GeV/c) 2 ] 

FIG. 10: The predicted Q 2 -dependence of meson cloud contributions to the helicity amplitudes 
A{. The dotted curves are from keeping only the 7r-loop in Eq. 0. The black dot corresponds to 
the absolute value of the bare helicities listed in Table. [H 



13 



V. Burkert and T.-S. H. Lee, Int. J. of Mod. Phys. E13, 1035 (2004). 

T.-S. H. Lee and L.C. Smith, J. Phys. G 34, 1 (2007). 

A. Matsuyama, T. Sato, and T.-S. H. Lee Phys. Rept. 439, 193 (2007) 

Feuster and U. Mosel, Phys. Rev. C 58, 457 (1998); Phys. Rev. C 59, 460-491 (1999); V. 
Shklyar, H. Lenske, U. Mosel, G. Penner, Phys. Rev. C 71, 055206 (2005). 
A. Usov and O. Scholten, Phys. Rev. C 72, 025205 (2005). 

A.V. Sarantsev, V.A. Nikonov, A. V. Anisovich, E. Klempt, and U. Thoma, Eur. Phys. J. A 
25, 441 (2005). 

G. Y. Chen, S. S. Kamalov, S. N. Yang, D. Drechsel and L. Tiator, Phys. Rev. C 76, 035206 

(2007) [arXiv:nucl-th/0703096] . 

S.N. Yang, J. of Phys. Gil, L205 (1985) 

S. Nozawa, B. Blankleider and T.-S. H. Lee, Nucl. Phys. A 513, 459 (1990). 
F. Gross and Y. Surya, Phys. Rev. C 47, 703 (1993). 

T. Sato and T.-S. H. Lee, Phys. Rev. C 54, 2660 (1996); Phys. Rev. C 63, 055201 (2001). 
S. S. Kamalov and S.N. Yang, Phys. Rev. Lett. 83, 4494 (1999). 

0. Krehl, C. Hanhart, S. Krewald, and J. Speth, Phys. Rev. C 60, 055206 (1999); C 62, 
025207 (2000). 

A. D. Lahiff and I.R. Afnan, Phys. Rev. C 66, 044001 (2002). 
M.G. Fuda and H. Alharbi, Phys. Rev. C 68, 064002 (2003). 
V. Pascalutsa and J. Tjon, Phys. Rev. C 61, 054003 (2000). 

W. T. Chiang, F. Tabakin, T.-S. H. Lee and B. Saghai, Phys. Lett. B 517, 101 (2001) 

B. Julia-Diaz, B. Saghai, T.-S. H. Lee and F. Tabakin, Phys. Rev. C 73, 055204 (2006). ' 

B. Julia-Diaz, T.-S. H Lee, T. Sato and L. C. Smith, Phys. Rev. C 75, 015205 (2007). 

R.A. Arndt, I.I. Strakovsky, R.L. Workman, Phys. Rev. C 53, 430 (1996); Int. J. Mod. Phys. 
A18, 449 (2003), Phys. Rev. C 66, 055213 (2002). 

D. Drechsel, O. Hanstein, S.S. Kamalov, and L. Tiator, Nucl. Phys. A645, 145 (1999); S.S. 
Kamalov, S. N. Yang, D. Drechsel, O. Hanstein, and L. Tiator, Phys, Rev, C 64, 032201(R) 
(2001). 

1. G. Aznauryan, Phys. Rev. C 68, 065204 (2003). 

S. Capstick and W. Roberts, Phys. Rev. D 49, 4570 (1994). 

P. Maris and C. D. Roberts, Int. J. Mod. Phys. E 12, 297 (2003). 

C. Alexandrou, G. Koutsou, T. Leontiou, J. W. Negele and A. Tsapalis, Phys. Rev. D 76, 
094511 (2007). 

D. Richards, larXiv:0711.2048l [nucl-th]. 

B. Julia-Diaz, T.-S. H. Lee, A. Matsuyama, T. Sato, Phys. Rev. C 76, 065201 (2007). 



CNS Data Analysis Center, George Washington University, http: / / gwdac.phys.gwu.edu] 

N. Suzuki, T. Sato, T.-S. H. Lee, MENU2007 proceedings, to be published. 

B. Julia-Diaz, T.-S. H. Lee, A. Matsuyama, T. Sato, L.C. Smith, in preparation. 

B. Julia-Diaz, D. O. Riska and F. Coester, Phys. Rev. C 69, 035212 (2004), [Erratum-ibid. C 

75, 069902 (2007)]. 



14 



